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Solving Legendre’s differential equation 


LEGENDRE's! differential equation has the form 
(1 —¢t?)y — 2ty + n(n + 1)y =0 (1) 


with the general solution 
y = aPp(t) + BQn(t), (2) 
where P,,(t) and Q,,(t) are linearly independent and are called LEGENDRE functions 
of the 1st- and 2nd-kind. P,,(t) are polynomials which is what we are interested in. 
Due to the factor 1—t? this equation can not be directly transformed into an analog 
computer setup as this would require a division, the denominator of which tends to 0. 
Instead equation (1) is rewritten as 


y —¢y —2ty+n(n+ly=0 
which is then solved for y, yielding 
yj = Py + Qty — n(n + 1)y. (3) 


This is insofar unusual as ¥ is not separated on the left side but also occurs on 
the right. In addition to that the terms 77% and 2tq would require three multipliers if 
implemented in a naive way. One multiplier can be saved by rewriting equation (3): 


y = t(ty + 2y) —n(n+ ly. 
To simplify notation, n* = n(n + 1) is introduced, finally yielding 
y = t(ty + 2y) —n*y. 


To select the desired solution P,,(t) for this DEQ suitable initial conditions are 
required. Table 1 lists the P,,(t) for 1 <n <5 with the resulting initial conditions for 
y(0) and y(0). Obviously the problem must be scaled at least with respect to y as the 
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Table 1: The first five P,,(t) solving LEGENDRE’s DEQ with the resulting initial conditions 


initial conditions nearly reach 2 for n = 5, thus giving rise to a scaling factor of 5 for 
y. 

n* must also be scaled — assuming a maximum value of n = 5 requires an additional 
scaling factor of a = 35 Some experimentation shows that i quickly approaches 
20 for n = 5 and t —> 1. Taking all these scaling factors into account yields the 
program shown in figure 1,2 the implementation of which on THE ANALOG THING 
is shown in figure 2. 

To generate P;(t) with this setup, the initial conditions as well as n* must be set 
according to table 1. It is (0) = 0 and y(0) = 2. With an overall scaling factor of 
539 the initial value must be set to 4(0) = sop © 0.047. n = 5 yields n* = 30 so 
that the potentiometer following the integrator yielding y must be set to 0.05. The 
operation time Is set to ns with the time scale factors shown. 

Figure 3 shows function y(t) = P(t) resulting from these values.? Experimenting 
with different initial conditions for y(0) and y(0) yields results of the general form of 
equation 2. 


1 ApRIEN-MARIE LEGENDRE, 18.09.1752-09.01.1833 
2The hat over the initial conditions denotes their scaled values. 
3This is an x, y-plot with 2 = t showing the flyback when THE ANALOG THING switched from OP to IC mode. 
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Figure 1: Analog computer program for solving the LEGENDRE DEQ 
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Figure 2: Actual setup of the problem on THE ANALOG THING 
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Figure 3: Ps(t) withO<t<1 
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